Plasmon polaritons in photonic superlattices containing a left-handed material 
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We analyze one-dimensional photonic superlattices which alternate layers of air and a left-handed 
material. We assume Drude-type dispersive responses for the dielectric permittivity and magnetic 
permeability of the left-handed material. Maxwell's equations and the transfer-matrix technique 
are used to derive the dispersion relation for the propagation of obliquely incident optical fields. 
The photonic dispersion indicates that the growth-direction component of the electric (or magnetic) 
field leads to the propagation of electric (or magnetic) plasmon polaritons, for either TE or TM 
configurations. Furthermore, we show that if the plasma frequency is chosen within the photonic 
(n(ui)) = zeroth-order bandgap, the coupling of light with plasmons weakens considerably. As light 
propagation is forbidden in that particular frequency region, the plasmon-polariton mode reduces 
to a pure plasmon mode. 
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Over the years, artificial complex materials have been 
increasingly used to shape and manipulate liglrbi^i 4 ^. 
The microstructuring of high quality optical materials 
yields remarkable flexibility in the fabrication of nanos- 
tructures, and allows for the tailoring of electromagnetic 
dispersions and mode structures to suit almost any need. 
Metamaterials 6 -^ ' 8 ' 9 ' 10 ' 11 ' 12 ' 13 , also known as left-handed 
materials (LHMs), are a remarkable example of such 
nanostructuring. Light propagation through metamatc- 
rials is characterized by a phase velocity opposite to the 
Poyinting vector, which corresponds to negative disper- 
sive electric and magnetic responses. 

The advent of metamaterials has opened a new era for 
optical devices, and has also given considerable thrust 
to the recent area of plasmonics, the study of plasmon 
polaritons. In metal-dielectric interfaces, for example, 
surface-plasmon polaritons are coupled modes that re- 
sult from resonant interactions between electromagnetic 
waves and mobile electrons at the surface of a metal or 
semiconductor. Such resonant surface-plasmon polari- 
tons may have much shorter wavelengths than that of 
the radiation, which enables them to propagate along 
nanoscale systems 1 -' 2 -' 3 -! 4 -, opening up a wide range of pos- 
sibilities for the construction of new optical devices, ft is 
well known, for example, that one of the most important 
features of surface-plasmon polaritons is to confine light 
to very small dimensions, yielding the merging of pho- 
tonics and electronics at the nanoscale. Furthermore, 
considering the dispersive character of the LHM's, to- 
gether with the enhanced optical magnetism exhibited 
by them, one might conjecture the existence of remark- 
able new phenomena such as the excitation of plasmon 
polaritons of a magnetic nature, that is, magnetic den- 
sity waves resulting from resonant interactions between 
the optical field and current densities at the metamate- 



rial. Ultimately, those developments could lead to in- 
creases in the resolution of microscopes, in the efficiency 
of LEDs, and in the sensitivity of chemical and biological 
devices 5 -^^. 

One-dimensional (ID) superlattices which al- 
ternate layers of positive and negative mate- 
rials have already exhibited many interesting 
properties 1 ^^ 1 ^^ ^! 2 ^ 3 - that are absent in 
superlattices composed solely of positive materials. 
In particular, the existence of a non-Bragg photonic 
bandgap, also known as a zeroth order gap, has been 
suggested^ 2 -, detecte d 13 ' 14 , and characterize d 15 ' 16 . In 
order to investigate the possibility of excitation of 
electric/magnetic plasmon polaritons, in the present 
work we study the oblique incidence of light on a model 
ID superlattice composed of layers A of air, and layers 
B of a doubly negative material. Layers A (width a) and 
B (width b) are distributed periodically so that d = a + b 
is the period of the superlattice nanostructure (cf. Fig. 

In the B layers, the electric and magnetic responses 
are dispersive and may assume negative values. If one 
neglects losses, they may be described by- 9 - ' 10 ' 11 ' 12 ' 13 



e B (u)) = E Q 



(1) 



where £_b(w) and hb(u) are the dielectric permittivity 
and magnetic permeability in slab B, respectively, one 
may choos o 10 ' 11 eo = 1.21 and [1q = 1.0, and the elec- 
tric/magnetic plasmon modes are at v — v e — 2 ^^; an< ^ 
which correspond to the solutions of 
0, respectively. 



v = v m = 



= and hb(u) 
We note that dispersions such as those in ^ hold 
in periodically LC loaded transmission lines 10 . These 
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negative-index systems were shown to exhibit good mi- 
crowave properties, with low loss and broad bandwidth^. 

We shall be interested in studying the properties of 
both the transverse-electric (TE: electric field parallel to 
the interface plane, see Fig. 1) and transverse-magnetic 
(TM: magnetic field parallel to the interface) polariza- 
tions of a monochromatic electromagnetic field of fre- 
quency u> propagating through a ID periodic system. 

In the case of a TE field, one may choose 



E(r, t) = E(z) exp [i (qx — ut)] e y 



(2) 



whereas, in the case of TM polarization, the magnetic 
field may be considered as 



H(r, t) = H(z) exp [i (qx — uot)] e y , 



(3) 



where we have assumed that the superlattice was grown 
along the z axis, q is the wavevector component along the 
x direction, and e y is the cartesian unitary vector along 
the y direction. 

Maxwell's equations lead to the following differential 
equations for the electric and magnetic amplitudes 
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where n(z) = \J z(z)yj \x(z) is the z-position dependent 
refraction index. 

In the sequel, ua and are the refraction indices, 
while fiA and /is are the magnetic permeabilities in lay- 
ers A and B, respectively. Note that, because of |T]), the 
riB refraction index in layer B is a function of the fre- 
quency, and may be a real positive, real negative or a 
pure imaginary number. Also, note that q, the incident 
wavevector component along the x direction, may be ob- 
tained as a function of the angle of incidence 9 = 9 a by 
q = ^n A s\n9. 

Equations ^ and may be solved by means of the 
transfer-matrix technique (see, for example, the work by 
Cavalcanti et al*£ and references therein). For ns a real 
number, and n 2 B —n 2 A sin 2 9 > 0, the procedure yields the 
TE polarization dispersion relation from the solution of 
the transcendental equation 



cos(kd) — cos(Q Aa) cos(Q sb) 

1 (Fa Fb 
F b F a 



sv&(Q A a) sin(Q B b). (6) 



In the formula, A; is the Bloch wavevector along the z di- 
rection which is the axis of the photonic crystal; fields in 
consecutive unit cells are related by the Bloch condition, 
i.e., by the phase factor e lkd . Qa and Qb are defined as 



}a = -n A \ cos6>| 
c 



-n A \ cos 6 A \, 



(7) 



s = s(z) fi = pt(z) 
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FIG. 1: (Color online) Pictorial view of the ID multilayer 
photonic superlattice with layers A and B in periodic arrange- 
ment, and the electric and magnetic fields for the TE-like and 
TM-like electromagnetic waves schematically shown. Note 
that, for normal incidence, the two polarizations are equiva- 
lent. 
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FIG. 2: (Color online) TE and TM dispersion relations v = 
v(k) in photonic periodic superlattices (y = ^-). Calculations 
were performed for air as slab A (ba = 1, A*a = 1), a = b = 12 
mm (also for a — 12 mm and b — 6 mm), and ui e /2n — 
LOm/2-K — 3 GHz for the Drude model [cf. Eq. (fTJ] in slab 
B. Dotted lines indicate the pure electric/magnetic plasmon 
modes at v — v e = „_" e ^- and v = u m = Solid and 



dashed lines correspond to 
angles, respectively. 



and 



7t/12 and 6 = incidence 



Qb = - y n% - n\ sin 2 9 = — \n B \ \ cos 9 B \ 



(8) 



where, F A = (Qa/^a), F b = {Qb/vb) and, in the 
last equality, we have made use of Snell's law ua sin 9 — 
us sin 9b- 

For us a real number, and n 2 B — n 2 A sin 2 9 < 0, the 
procedure yields the TE polarization dispersion relation 
from the solution of the transcendental equation 



cos(kd) = cos(Q ao.) cosh(Q sb) 

2\F b F a 
where Qa is still given by |J7J, but 
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sm(Q A a) sinh(Q s &). (9) 
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FIG. 3: (Color online) TE dispersion relations v = v(k) in 
photonic periodic superlattices (y — ^-) as functions of vari- 
ous incidence angles. Calculations were performed using the 
same parameters as in Fig. 2. Dotted lines indicate the 
pure magnetic plasmon modes at v = v m = u *2— , dashed 
lines correspond to 6 = 0, and solid lines correspond to (a) 
9 = tt/12, (b) 9 = tt/6 and (c) 9 = tt/3 incidence angles, 
respectively. 
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FIG. 4: (Color online) TE and TM dispersions in a photonic 
periodic superlattice with a — b = 12 mm. Numerical calcu- 
lations were carried out for air as slab A (ea = 1 and ha = 1). 
For slab B, es and {ib were chosen by taking plasma frequen- 
cies in the Drude model as uj e /2ir = 1.5 GHz and oj m /2n = 2.5 
GHz, respectively, in (a) and (b), whereas in (c) and (d) we 
have taken u e /2n = 2.5 GHz and uj m /2iT — 1.5 GHz. Note 
that, in panel (b), the electric plasmon frequency v e = 2-^/6^ 
is in the {n(u)) = zeroth-order bandgap region, leading 
to a basically flat 9 — tt/12 plasmon-polariton band with a 
frequency essentially equal to v e for the TM dispersion. A 
similar situation is observed for the TE dispersion in (c), as 
v m = 2 ^L_ falls in the zeroth-order gap region. Solid and 
dashed lines correspond to 6 = 7r/12 and 9 = incidence 
angles, respectively. 



Moreover, if n 2 B < 0, © is still valid, with Qa still given 
by ©, but 

Qb = ^\Jn 2 BI -n 2 A sin 2 6>, (11) 

where n bi is the imaginary part of ub ■ 

For the TM polarization, the transcendental equations 
© and ©, as well as the definitions of Qa and Qb for 
the different cases are also valid, provided one replaces \ia 
by ca and fiB by cb, where €a and es are the dielectric 
permissivities in layers A and B, respectively. 

Figure [21 shows the calculated dispersions (v = in 
the cases of TE [Figs. Ha) and He)] and TM [Figs.Hb) 
andHd)] polarizations, for different layer widths. Note 
that v e and v m fall outside the {n(uj)) = zeroth or- 
der bandgap. For 6 — (dashed curves), the second 
photonic band, just above the (n(uj)) = zeroth or- 
der gap, is a pure photonic mode and, for 9 — tt/12, 
there is a coupling between the radiation field and a plas- 
mon mode which leads to a pair of plasmon-polariton 
modes. In other words, while for normal incidence the 
band just above the zeroth order gap is a pure photonic 
mode, for oblique incidence (9 ^ 0) the electromagnetic 
field + plasmon interaction leads to a pair of coupled 



modes. We emphasize that, for 9 = 0, the electric- and 
magnetic-field components along the z growth direction 
are null and there is no longitudinal wave propagation. 
On the other hand, for 9 ^ 0, in the case of the TM 
(TE) configuration, the electric (magnetic) field has a 
component on the z direction, which results in the ex- 
citation of longitudinal electric (magnetic) waves. For 
9 = tt/12, Figure [J] illustrates that the pair of plasmon- 
polariton bands are asymptotic to the v = v e and v = v m 
pure electric/magnetic plasmon values. It is clear from 
Fig. [2] (a), for example, that at small values of k, the 
lowest 9 = tt/12 plasmon-polariton mode behaves like 
an electromagnetic photonic wave. As k increases, the 
dispersion bends over and reaches the limiting value of 
the pure magnetic plasmon. On the other hand, the sec- 
ond 9 = tt/12 plasmon-polariton branch behaves like a 
magnetic plasmon wave, at small k, and like an electro- 
magnetic photonic mode as k increases. The same com- 
ments are valid for the other panels in Fig. [5] [of course, 
in the cases of Figs. Hb) andH°0> the electric plasmon 
mode is the one involved] . By comparing the results for 
normal and oblique incidence, it is clear that, for 9 =/= 0, 
resonant plasmon-polariton waves are excited by the cou- 
pling of the electric (or magnetic) plasmon modes with 
the incident eletromagnetic field. One may infer that 
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the plasmon-polariton waves are excited by the magnetic 
field, in the TE case, or driven by the electric field, in 
the TM case. This is consistent with the fact that those 
longitudinal waves, in the TE case, are asymptotic to the 
v m pure magnetic plasmon frequency whereas, in the TM 
case, the plasmon-polariton bands are asymptotic to the 
v e pure electric plasmon mode. 

The calculated dispersions, for TE modes, as a function 
of the incidence angles are displayed in Fig. [3] We no- 
tice that the essential features persist irrespective of the 
incidence angle, although the Bragg gaps widens consid- 
erably for higher incidence angles. Also, Fig. Q] clearly in- 
dicates that, by choosing the resonant plasmon frequency 
within the (n(u>)) = zeroth order bandgap, the coupling 
of light to plasmon modes, for 8^0, essentially disap- 
pears leading to a basically pure (electric or magnetic) 
plasmon mode. This is a consequence of the fact that 
the energy of the incident electromagnetic wave lies in a 
forbidden energy gap region and, therefore, the coupling 
of the incident light with plasmons is expected to be quite 
weak. 

Summing up, we have studied light propagation 
through a ID superlattice composed of alternate lay- 
ers of a positive constant material and a negative dis- 
persive material, and verified the appearance of cou- 
pled plasmon-polariton modes of electric and magnetic 
natures. The photonic dispersion was calculated from 
Maxwell's equations using the transfer-matrix approach, 
with the LHM modeled by Drude-like dielectric and mag- 
netic responses. In the case of oblique 9^0 inci- 
dence, the photonic dispersion indicates that the mag- 
netic or electric field in the frequency region around the 



(n(u))) = zeroth order bandgap leads to coupled mag- 
netic or electric plasmon-polariton modes for the TE 
and TM configurations, respectively. Moreover, present 
results show that the coupling of light with plasmons 
is weakened by choosing the plasma frequency inside 
the zeroth order gap. As light propagation is forbid- 
den in that particular gap-frequency region, the coupled 
plasmon-polariton mode becomes essentially a pure plas- 
mon mode. This feature permits one to select which type 
of plasmon polariton (i.e., electric or magnetic) one is 
willing to excite, by choosing the magnetic or electric 
plasmon frequency within the zeroth order gap. Here, 
one might conjecture that the possibility of excitation of 
magnetic plasmon polaritons might provide a new avenue 
to the implementation of novel techniques and devices 
based on the interplay of photonics and magnetism at 
the nanoscale, similarly to the interplay of photonics and 
electronics. Therefore, the implementation of plasmonic 
chips for high data rate processing or efficient sensing ap- 
plications might also be obtained based on the excitation 
of magnetic plasmon polaritons. Finally we have studied 
an ideal system in which losses have been neglected, and 
future work including such effects is certainly forthcom- 
ing. 
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